Electronic structure of V2O3 : Wannier orbitals from LDA-iVMTO calculations 
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Using muffin-tin orbital (MTO) based iVMTO-downfolding procedure within the framework of 
local density approximation, we construct the Wannier orbitals for the ti g manifold of bands in 
V2O3 in the paramagnetic phase. The real space representation of the one-electron Hamiltonian 
in the constructed Wannier function basis shows that, contrary to the popular belief, the in-plane 
hopping interactions are as important as the vertical pair hopping. Following the language of Di 
Matteo et.al. [Phys. Rev. B 65 054413 (2002)], this implies, the problem of V2O3 falls in the atomic 
regime rather than in the molecular regime. We have also repeated our construction procedure in 
the low temperature monoclinic phase, for which the changes in hopping interactions are found not 
to be dramatic. 

PACS numbers: 71.20.-b, 71.15.Ap, 71.20.Bc 
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I. INTRODUCTION 

V2O3 has been in focus of attention since 1969 when its 
unusual phase diagram was discovered 1 -. At low temper- 
ature, pure V2O3 is an antiferromagnetic insulator (AFI) 
with a monoclinic, slightly distorted corundum struc- 
ture, a complicated magnetic order, a moment of 1.2 \xb-, 
and a gap of 0.66 eV 2 . At T/v=154K it transforms to 
a corundum-structured, paramagnetic metal. Upon sub- 
stituting V by Ti or by application of pressure, the Neel 
temperature decreases and the antiferromagnetic phase 
vanishes above 5% Ti. Substitution of V by Cr, on the 
other hand, causes the Neel temperature to increase and 
reach 180 K for 1.8% Cr. For higher concentrations, Tjv 
stays constant and the transition is to a paramagnetic 
insulator. For Cr concentrations between 0.5 and 1.8% 
there is a second phase transition, which upon increas- 
ing temperature, or Cr concentration, is from a para- 
magnetic metal (PM) to a paramagnetic insulator (PI). 
This transition is isostructural, ends at a critical point, 
(T c , y c ) = (400 K, 0.5% Cr) , and has been considered the 
classic example of a Mott-Hubbard transition. 

In the high-temperature corundum structure (FIG. 
[1]), all vanadium ions are equivalent and surrounded by 
nearly perfect oxygen octahedra. Since the covalent O- 
V pda interaction is stronger than the pdir interaction, 
the more antibonding Vd-like e g level 3 - lies above the 
less antibonding V rf-like ti g level and, as a consequence, 
the electronic configuration of V2O3 is Vif ff - Now, the 
three-fold degenerate ti g level is split into an upper a\ g 
and a lower, doubly degenerate e g level by a trigonal dis- 
tortion, which mainly consists of a slight displacement 
of the vanadium ions along the vertical three-fold axis, 
away from the centers of their octahedra, so that the dis- 
tance between a vertical vanadium pair (V 4-1 or 2-5 in 
FIG. [T|) is slightly longer than the distance between the 
centers of the two octahedra. 



There have been many attempts to explain these 
metal-insulator transitions and the spin structure of the 
antiferromagnetic insulating phase. The careful analysis 
presented in 1978 by Castellani, Natoli, and Ranningcr 
(CNR) 4 resulted in a model which remained undisputed 
for over twenty years: Since the a\ g orbitals have d,3 Z ^-i 
character and point towards each other, yielding a strong 
dda-like hopping integral, these orbitals on each verti- 
cal pair form bonding and antibonding levels which are 
split by more than twice the a\ g -e g crystal-field split- 
ting and by more than the on-site Coulomb interaction. 
As a consequence, one electron per vanadium is used 
to form a spin-singlet, chemical bond between a verti- 
cal pair. The other electron enters the doubly degener- 
ate, localized orbitals. Since the integrals for hopping 
from- and between e g orbitals are relatively small, the 
on-site Coulomb repulsion leads to an 5=1/2 state and 
may order the occupied orbitals in a way consistent 
with the observed spin structure of the low-temperature 
antiferromagnetic insulator. In this structure, the spins 
on vanadium pairs in the z and x directions are aligned 
ferromagnetically, and those on pairs in the other two di- 
rections are aligned antiferromagnetically. This requires 
an orbital order in which an integral for hopping in the 
x and z directions between occupied orbitals is con- 
siderably smaller than between occupied and unoccupied 
Cg orbitals. This was reviewed and discussed by Rice 5 
using Kugel-Khomskii's general description 6 of the cou- 
pling between orbital and spin degrees of freedom. The 
CNR model also led to a half-filled, one-band (the low- 
est e^-band) Hubbard Hamiltonian to serve as the sim- 
plest possible electronic model for V2O3. This model was 
solved by Rozenberg et. al£ using the dynamical mean- 
field approximation (DMFT) 8 and found to describe the 
metal-insulator transition. However, the polarized x-ray 
absorption experiment of Park et. al£, corroborated with 
multiplet calculations showed that V 3c? 2 ions are in the 
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high spin (5 = 1) state rather than in 5 = 1/2 state 
and the orbital occupation, which is different in different 
phases, is an admixture of e g e g states with e g a\ g config- 
urations. From that, they concluded that neither V2O3 
problem can be mapped onto a single-band Hubbard 
model, nor the projecting out of a\ g orbitals by means 
of molecular orbital formation as was done by CNR^ is 
justified. Ezhov et. aVs calculations 10 within the local 
density approximation (LDA)+U scheme showed the im- 
portance of Hund's rule exchange giving rise to 5 = 1 
model but with only e g occupancy implying no orbital 
ordering. However, they succeeded in correctly predict- 
ing the low-temperature magnetic structure which was 
attributed to be stabilized by the monoclinic distortion. 
Nevertheless, the issue associated with the orbital order- 
ing remained which apparently showed up its presence in 
several different experimental observations^. To recon- 
cile the 5 = 1 and the orbital ordering aspect, Mila et. 
alH and Di Matteo et. alM- subsequently proposed two 
different correlated model of c — axis pair states incorpo- 
rating dynamical mixing of e g e g and e g ai g states with 
5 = 1 spin configuration on each of the sites. In recent 
years, combined with the LDA, DMFT calculations have 
been carried out«. 

While, it is now generally accepted, a realistic theory 
of V2O3 must take into account the complicated elec- 
tronic structure of the system, there has been no seri- 
ous attempt to realistic modeling of the electronic struc- 
ture of V2O3 since the early work of CNRi which was 
crude in its various approximations and was partly semi- 
empirical. The starting point of several of the many-body 
model-based calculations seem to be the vertical V-V pair 
model which is considered to be the predominant build- 
ing blocks — the validity of such assumptions need to 
be re-examined in the context of accurate tight-binding 
(TB) modeling of V 2 3 . 

In recent years MTO based N—th order MTO method, 
namely ATVITO metho d 15 i 16 has been introduced and im- 
plemented. The method goes beyond the scope of the 
standard linear MTO (LMTO) method, in defining an 
energetically accurate basis set with a consistent descrip- 
tion throughout the space of MT spheres and the inter- 
stitial. An important feature of the 7VMTO method is 
the so-called downfolding technique which provides an 
useful way to derive few-orbital Hamiltonians starting 
from complicated full LDA Hamiltonian by integrating 
out degrees of freedom not-of- interest. This procedure 
naturally takes into account the renormalization effect 
due of the integrated-out orbitals by defining energy- 
selective, effective orbitals which serve as the Wannier 
or Wannier-likc orbitals for the few-orbital Hamilto- 
nian in downfolded representation. The method pro- 
vides a first-principles way of deriving the single-particle 
model Hamiltonian and direct generation of Wannier 
functions without any fitting procedure giving rise to 
an unique scheme that has the deterministic nature of 
first-principles calculations added to the simplicity of 
model Hamiltonian approaches. The method has proved 



to be extremely successful in deriving model Hamilto- 
nians for systems such as high-Tc cupratesil, double 
perovskites^, low-dimensional quantum spin systems^. 
This approach of direct generation of Wannier functions 
may be contrasted to that of construction of Wannier 
functions out of the calculated Bloch functions. Re- 
cently Anisimov et at implemented a similar method for 
use in LDA+DMFT= and Solovyev proposed a general 
LMTO-based procedure for constructing effective lattice 
fermion models 2 —. With other local-orbital basis sets, 
somewhat similar techniques can be used^ 2 -, but in case 
not all basis functions are well localized, e.g. for the set 
of bare LMTOs^, the Wannier functions obtained for 
the correlated bands may not be sufficently localized for 
the corresponding on-site-t/-Hamiltonian to be realistic. 
For those cases, more complicated procedures for obtain- 
ing for instance those Wannier functions which minimize 

the spread^ ^|r — (r)| 2 ^) , or those which maximize the 
Coulomb self-energy, has been used.— ^ 

In this paper, we aim to provide an accurate 
tight-binding description of V2O3 by constructing the 
Wannier-like functions for the ti g Hamiltonian employing 
the ATVITO methodology. In this context, this methodol- 
ogy has recently been applied 2 ^ for tight-binding model- 
ing of aig bands of ferromagnetic LDA+U calculations to 
investigate the role of vertical pair from band-structure 
point of view. In those specific calculations, primary in- 
terest was to estimate the c — axis intra-pair hopping 
matrix element compared to inter-pair hopping matrix 
elements for a\ g bands and therefore, LDA+U rather 
than LDA was chosen as the basis of calculations, which 
provides nice separation of a\ g and bands. However, 
from point of view of input to realistic many-body cal- 
culations it is more suitable and preferable to start with 
LDA-derived Hamiltonians and the full t2 g Hamiltonian 
since the a\ g and states are both equally important as 
seen in experiment and the hybridization effect between 
a\ g and e g should be taken in account. In the present 
work, we have therefore chosen LDA as the basis of our 
calculations. The present paper in that respect, should 
be considered as a more detail paper for the tight-binding 
modeling of V2O3 system. The tight-binding parame- 
ters derived in this paper in a rigorous, first-principles 
manner will be useful as an input to many-body vari- 
ational calculation like that of Di Matteo et. al^ as 
will be discussed in section III D. The parameters can 
also be used for the many-body LDA+DMFT calcula- 
tions. The LDA+DMFT calculations report in Rcf. 28 , 
has been carried out using the AMTO Wannier func- 
tion implementation of LDA, presented in this paper. 
These calculations 2 ^ showed the importance of correla- 
tion assisted dehybridization of a\ g and e g in the de- 
scription of the correlated electronic structure of V2O3 
and its metal-insulator transition. Calculations within 
such LDA+DMFT framework has been also used to ex- 
plore the comparison of doping, temperature and pres- 
sure route to metal-insulator transition in V 2 0^ and 
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to study multi-orbital effects in optical properties of 

In the following, in section II we discuss the crystal 
structure and orbital symmetry aspects in V2O3. Sec- 
tion III involves description of the results and discus- 
sion. This section is divided into several sub-sections. 
In the sub-section A we present the high-energy part 
of the LDA band-structure of rhombohedral, undoped 
V2O3 in ambient pressure. In the sub-section B we ex- 
plain how downfolding within the AMTO method can 
be used to construct truly minimal basis sets which pick 
out selectively O-p, V-t2g, Ve 9 or V-s bands. We also 
present the TVMTOs - the members of such truly mini- 
mal sets of V-t2g and Ve 9 . With this tool at hand, in 
sub-section C we zoom in on the LDA t 2g bands, where 
due to lowering of symmetry induced by trigonal distor- 
tion we switch on from the t 2g representation to more 
appropriate a\ g — e g representation. We present the min- 
imal set constructed out of a\ g and e g 's, their Wannier 
functions, the hopping integrals and the comparison with 
existing results. In sub-section D, we discuss the validity 
of vertical pair model and the molecular orbital based 
approaches in light of 7VMTO derived hopping integrals. 
We restricted our study to high temperature paramag- 
netic phase of pure V2O3 in rhombohedral, corundum 
structure until the sub-section E, where for the sake of 
completeness of our study, we also discuss the hopping 
integrals in the Cr-doped V2O3 (sub-section E) and in 
the low-temperature monoclinic structure (sub-section 
F) . Finally, we conclude in section IV with summary and 
outlook. The essential details about the ATVITO method, 
which will be used for construction of localized Wannier 
orbitals and truly minimal, downfolded basis sets can be 
found in the Appendix A. 



structure. This structure consists of hexagonal packing 
of the oxygen atoms, and the vanadium atoms occupy- 
ing 2/3 of the octahedral cation sites. The basic features 
of the corundum structure is shown in FIG. [T] The im- 
mediate surrounding of the V atoms provided by oxygen 
atoms has approximate octahedral symmetry. The VOg 
octahedra face share along the vertical direction forming 
V-V vertical bonds, while they edge share forming lay- 
ers of honeycomb lattice, giving rise to three-dimensional 
network with overall rhombohedral symmetry. In the re- 
sulting corundum structure, which has R3c space group 
symmetry, each primitive rhombohedral unit cell con- 
tains two V2O3 formula units while the non-primitive 
hexagonal unit cell contains six V2O3 formula units. The 
experimentally determined structure^ of ambient pres- 
sure, pure V2O3, in the hexagonal setting, with lattice 
constants a#=4.952 A and ch =14.003 A, and V and O 
atoms occupying the Wyckoff positions (12c) and (18e) 
with internal parameters zy — 0.34630 and xo — 0.31164, 
yield V-0 bond lengths in the range 1.971 A- 2.049 A. 
The nearest-neighbor V-V distances within the hexago- 
nal layers are 2.882 A while that along the vertical di- 
rection is about 6 % shorter (2.697 A). As is evident 
from the structural figure as well as from the internal 
parameter value, the V atoms are displaced from their 
ideal positions where V atoms in the hexagonal layers 
would have been co-planer. This displacement causes V 
atoms to move away from the center of the octehedra, 
giving rise to three long and three short V-0 bonds. The 
arrangement of V atoms along the hexagonal z-axis can 
be derived from an ideal chain structure by introducing 
vacancies at every third site. 



II. CRYSTAL STRUCTURE AND SYMMETRY: 



As mentioned in the introduction, the high temper- 
ature paramagnetic phase of V2O3 has the corundum 



The results reported in the following, for the high- 
temperature structure, are all carried out with rhombo- 
hedral unit cell defined by the primitive lattice transla- 
tions in a Cartesian system with z-axis pointing along the 
vertical V-V bond and the x-axis chosen as the projection 
of V1-V2 onto the xy plane: 
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where c/a = 1.633 and a = 2.859 A. The 3-fold axis 
pointing along the z— direction are given by Ti + T2 + 
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FIG. 1: (Color online) The crystal structure of V2O3 in the high-temperature paramagnetic phase. The larger (green) circles 
indicate the V atoms. The smaller (violet) circles surrounding the V atoms are oxygens, showing the octahedral co-ordination. 
The VO6 octahedra face-share along the vertical direction while they edge-share within the hexagonal layers. The unprimed 
and primed co-ordinate systems represent the rhombohedral co-ordinate system (z-axis pointing along the vertical V-V bond 
and x-axis chosen as the projection of VI- V2 in the xy plane) and the oxygen-based octahedral co-ordinate system (z -axis 
pointing along 06-05, x'-axis pointing along 03-01) respectively. 



The high-symmetry points on the Brillouin zone (BZ) 
as shown in FIG. [2] are given by, 

2 points on the 3-fold axis and at the center of 
regular hexagonal faces of the BZ : 



1 2tt 
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6 points at the center of hexagonal faces of BZ with 2 
short and 4 long edges: 
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and 6 points at the center of rectangular faces: 
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The perfect octahedral crystal field surrounding of the 
V ions split the 3d energy levels into two-fold degenerate 
e g levels, d 3z ,2_ 1 , & x ,i_^,i and three- fold degenerate i 2 g 
levels, d x ' y ', <i x ' z ', d„' z ', where the primed co-ordinate 



■5 




FIG. 2: The rhomb ohedral BZ showing the high-symmetry 
points. 

system refers to the octahedral co-ordinate system with 
x ,y and z pointing along 03-01, 04-02 and 06-05 
[see FIG. [TJ. However, the oxygen octahedra surround- 
ing of the V ion in V2O3 is not quite perfect, but has the 
trigonal distortion. This trigonal distortion of the octa- 
hedral environment of the V site and influence of non- 
cubic arrangement of more distant V ions in the lattice 
lowers the symmetry from octahedral Oh group to 
group, resulting into further splitting of the t 2g complex 
into singly degenerate a\ g and two- fold degenerate e^. 
Starting from three congruent t2 g orbitals, d,^' , d x < z <, 
d y i z > which can be derived from each other by a counter- 
clockwise rotation of 27r/3 around the three- fold z-axis, 
the dig and orbitals are generated as: 



d m = ±(d x , z ,+d x , y ,e 2 ^ 3 + d y , z/ e-^ m / 3 ) 



dig : 771 = 0, : m = ±1 
e* 1 - V2Im d 1 = l/V2(d xW -d y , z ,) 
e*,2 = V2Re d 1 = v^T&W " W^y + d y , Z Q) 



Further, transforming to rhombohcdral, unprimed co- 
ordinate system x, y, z, the three orbitals transforming as 
the dig and representations are given by^, 



aig ■ d 3z 2_ 1 

e* 1 = y/2/3d xy + 1/V3d xz 

e^, 2 = -^/~2/3d x 2_ y 2 - l/V3d yz 
The high-temperature Cr-doped V2O3, which is para- 
magnetic, insulator in nature, retains the corundum crys- 
tal structure, with lattice constant expanding to 4.998 
A and c/a dropping to 2.78 for substitution of approx- 
imately 1 % Cr, compared to undoped, metallic phase 
discussed above. This causes expansion of the vertical 
and all nearest-neighbor basal V-V bonds by 1.8-14 %. 
Upon doping with Cr, the Wyckoff positions of V and O 
also change. For (Vo.962Cro. 038)203 they become zy — 
0.34870 and x = 0.30745, yielding V-0 bond lengths in 
the range 1.976 A- 2.061 A. 

In the low-temperature AFI phase, the crystal struc- 
ture is further distorted from the corundum structure to 
monoclinic. This distortion causes the tilting of the ver- 
tical V-V bond by 1.8° towards the positive side of the 
x— axis and breaks the three-fold rotational symmetry, 
resulting into monoclinic crystal of symmetry I2 /a with 33 
a m = 7.255 A, b m = 5.002 A, c m = 5.548 A, p = 96.752° 
and four formula weight per unit cell. As a consequence, 
the vertical V-V bond length increases slightly from 2.697 
A to 2.745 A, one of the V-V bond (V1-V2) within the 
hexagonal layer elongates to 2.986 from undistorted bond 
length of 2.882 A, while the other two (VI- V3 and VI- 
V3 ) remain essentially same with bond lengths 2.862 A 
and 2.876 A. In the low-temperature magnetic struc- 
ture which is rather unusual, the vertical V1-V4 bond 
and the VI- V2 bond along the x axis becomes ferro- 
magnetic while the other two basal bonds become an- 
tiferromagnetic. Upon distortion, the oxygen octahedra 
also becomes slightly skewed about the central V atom, 
while the average V-0 bond length remains practically 
unaltered. In FIG.[3]we show the low-temperature mon- 
oclinic structure together with high-temperature corun- 
dum structure. The low symmetry crystal field in the 
monoclinic phase, further lifts the degeneracy between 
two e^ orbitals and mixes the a\ g and orbitals on the 
same site. 



III. RESULTS AND DISCUSSIONS about 16 eV around the Fermi level (set as zero in the fig- 

ure). The bands are plotted along the various symmetry 
A. Corundum V2O3 : LDA band-structure directions of the rhombohedral BZ, shown in FIG. [H The 

results are obtained with self-consistent potentials gener- 

FIG. U] shows the LDA one-electron band-structure of 
corundum V2O3 in PM phase, over an energy range of 



FIG. 3: (Color online) The crystal structure of V2O3 in the low-temperature monoclinic phase. For comparison, the high- 
temperature corundum structure is also shown in the same figure. The structures have been chosen to match at the central 
point of the vertical VI- V4 bond. As in FIG. [T] the larger (red for corundum and green for monoclinic) circles indicate the V 
atoms. The smaller (blue for corundum and violet for monoclinic) circles surrounding the V atoms are oxygens. 



ated out of the tight-binding LMTO calculation within 
the atomic sphere approximation (ASA}21. The details 
of the computation may be found in Appendix B. von 
Barth and Hcdin parametrization-^ has been used for the 
LDA exchange-correlation potential. The band-structure 
results presented in the FIG. 0] are obtained with stan- 
dard set of nearly-orthonormal LMTO's, whose accuracy 
is good enough for describing the high-energy features of 
the band-structure. The bands are in good agreement 
with the linear-augmented-plane-wave (LAPW) result of 
Mattheiss2M2. 

Plotted bands are the orbital-projected bands or the 
so-called fat-band in the sense that the fatness of the 
bands in each panel is the weight of the indicated orbital 
in the wave- function. The co-ordinate system is chosen as 
that of the oxygen-based octahedral co-ordinate system 
with the z -axis pointing along the 06-05 direction and 
x-axis pointing along the 03-01 direction. As is seen, 
the low-lying bands below -3 eV is predominantly of oxy- 
gen character. With the choice of octahedral co-ordinate 
system, the V 3d splits into ti g and e g manifolds. 12 t2 g - 
like bands (since there are 4 V atoms in the unit cell of 
the primitive rhombohedral unit cell with 3 t^g orbitals 
on each V ion) lying lower in energy compared to e g -like 



bands cross the Fermi level, spanning an energy window 
from about -1.5 eV to 1.5 eV. The crystal field split eg- 
like bands lye high up in energy from about 1.7 eV to 
4 cV separated from the t2 g manifold by a small energy 
gap of about 0.2 eV. V-s dominated states lye further 
high up in energy starting from about 4.5 eV. 

In the oxygen-projected band-structure, we notice, in 
addition to predominant fatness associated with oxygen- 
dominated bands lying below -3eV, the fatness associ- 
ated also with V-d dominated bands. Similarly in V-d 
projected band-structure we notice the presence of char- 
acter in 0-p bands, which is born in by the V-d — 0-p 
hybridization. It is this V-d — 0-p hybridization, that 
moves the V-d and 0-p dominated states far apart from 
each other with oxygen bands fully occupied and V-d 
bands mostly empty. Due to the different orientation, 
e g orbitals hybridize more strongly with 0-p forming di- 
rected pda bonds while the <2g orbitals bond less strongly 
giving rise to pdir bonds. This is evident from the fat- 
band plots which shows that the upper part of the V-d 
bands - the energy region dominated by the V-e g and the 
lower part of the 0-p bands has the most mixing in. We 
also notice the significant hybridization between V-s and 
0-p degrees of freedom. 




FIG. 4: LDA band-structure of V2O3 in the high-temperature corundum structure, plotted along the symmetry directions of 
the rhombohedral BZ. The BZ is shown in FIG. [2] The fatness associated with each band is proportional to the character of 
the orbital indicated at the top of each panel. Zero of the energy is set at the Fermi level. 



B. Corundum V2O3: downfolded few-orbital 
band-structure and Wannier-functions 

In the following, we demonstrate the application of 
NNLTO-downf aiding technique to produce truly minimal 



basis sets which may be chosen to span selected bands 
with as few basis orbitals as there are bands. This is 
illustrated by constructing truly minimal basis sets for 
O-p, V-i 2g , V-e g and Vs. 



In the first three panels of FIG. we show the bands 
obtained by using truly minimal sets, either O-p or V-i 2 g 
or V-e g in solid lines as compared to full LDA band- 
structure in dotted lines. The basis sets for the band- 
structure calculations shown in solid lines, which we call 
as downfolded bands, contain as many orbitals as the 
number of bands - hence is the name truly minimal ba- 
sis set. The NMTO-downfolding procedure enables one 
to construct a set of O-p or V-i 2g or V-e g muffin-tin or- 
bitals of order N, AMTO, which span the O-p-like or 
V-i 2g -like or V-e g -like bands - and no other bands - with 
arbitrary accuracy as N increases. Such a set is exact 
for the energies, eq, —• ,ew, chosen for its construction. 
As is seen in FIG. [5j three energy points were used for 
the construction of O-p, V-i 2g and V-e g minimal sets, so 
the MTOs are of order N = 2, i.e. they are quadratic 
MTOs or QMTOs. Since the 18 O-p-like, 12 i 2g -like and 
8 e g -like bands are isolated from the above and below- 
lying bands, the AMTO set obtained by making the en- 



ergy mesh, eo,....,ejv, finer and finer will converge to 
the Hilbert space spanned by any set of Wannier func- 
tions. In other words, the symmetrically orthogonalizcd 
set of converged AMTOs is a set of Wannier functions. 
As is seen, already with choice of three energy points, 
the downfolded bands are indistinguishable from the full 
LDA bands in the region of interest spanned by O-p, V- 
i 2g and V-e g respectively - so, they are converged in the 
above-mentioned sense. The corresponding AMTOs are 
localized by construction as explained in the Appendix, 
but they are not quite orthogonal. These truly minimal 
AMTO sets, therefore must be symmetrically orthogo- 
nalized in order to become a set of localized Wannier 
functions. In FIG. [6] we show one of the three congru- 
ent orbitals of such a f 2g AMTO set, namely d x > ' and 
two orbitals of the e g set before orthogonalization. Only 
the central part of the orbitals have d x ' y i, d x ,2_ y ,2 or 
d 3z '2_ 1 character. In order to describe the hybridization 
with the O-p and the hybridization between V-t2 g and 
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V-t2g 



V-eg 



V-s 





FIG. 5: LDA band structure of corundum-structured V2O3 in various energy ranges. The solid lines in different panels show 
the bands obtained with the truly minimal (downfolded) 0-2p, V-i22 9 , V-e 9 and V-4s iVMTO basis in comparison to those 
obtained with full NMTO basis (shown in dotted lines). Apart from the high-energy range in the last panel, where only V-4s 
ATVITO's are used to form the truly minimal set, the downfolded bands in various panels are indistinguishable from the bands 
in the full ATV1TO basis, within the respective energy range of interest. The NMTO energy points, e„-s, spanning the region of 
interest are shown on the right-hand side in each panel. 



V-e 9 within the V-d manifold, the O-p character and, 
V-e 9 character for the case of d x ' 1 and V-t2 g character 
for the case of d 12 « or i n ,, are folded into the 
tails. In fact, all the partial wave characters other than 
the respective active characters, d x i y > in the left panel, 
d 12 12 in the middle panel or cL ,2 in the last panel, 
are folded down in the tails. We see the strong pda anti- 
bonds in the plots of e g AMTOs and relatively weak pdir 
anti-bonds in the plot of t 2g NMTO. 

The last panel in FIG. [5j deals with the more difficult 
case, where the chosen bands, namely the V-s bands over- 
lap with the other high-lying bands, e.g. V-p, O-d bands. 
As is seen in the figure, even for such a difficult case of 
bands of interest overlapping with other bands, it is pos- 
sible to pick out the selected bands - in the present case 
four V-s bands arising from four V atoms in the unit cell. 
With chosen two energy points, the downfolded bands 
(shown in solid lines) differ from the full LDA band- 
structure (shown in dotted lines). Nevertheless, we sec 
that the bottom part of the V-s derived bands has been 
reproduced quite well over an energy range of about 2 
eV with merely two energy points. As expected, increas- 



ing the number of energy points improves the agreement. 
The Wannier functions of such a complex of bands which 
is overlapping with other band complexes are ill-defined - 
the corresponding set of orthogonalized AMTOs, is there- 
fore, the set of Wannier-like functions. 

As an illustrative purpose for the LMTO practitioners 
and to appreciate the improvements within the AMTO 
procedure, in FIG. [7] we also show the downfolded O-p, 
V-t2 g , V-e 9 and V-s computed within the framework of 
LMTO, using the standard TB-LMTO code where the 
nearly orthogonal LMTO's are used for producing the 
truly minimal basis sets. As shown in the figure, the 
method works to a certain level of accuracy, provided 
all the I- and independent e„'s are put in the energy 
region of interest and not at the center of gravity of the 
respective occupied manifold as is done during the self- 
consistent loops of LMTO. However, within the LMTO 
scheme, the accuracy of the downfolded bands compared 
to the full basis band structure is not up to the level of 
satisfaction and importantly the generated minimal basis 
does not have the desired Wannier-like description. 



J 



Corundum V2O3: Tight-Binding Hamiltonian 
corresponding to t2 9 bands 



t2 g manifold that gets partially filled with two V elec- 



A reasonable approach in the tackling the V?Ot prob- 
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FIG. 6: (Color online) One of the three congruent orbitals 
of the truly minimal t2 g NMTO set and the two orbitals of 
the truly minimal e g NMTO set corresponding to downfolded 
bands in FIG. [5] Shown are the orbital shapes (constant- 
amplitude surfaces) with the lobes of opposite signs labeled 
by red and blue respectively. iVMTO's are localized by con- 
struction: An orbital of e.g. t2 g set is confined by the condi- 
tion that it has no t2 g character on other V atoms but may 
have O-p, V-e 3 characters, as is evident from the plots. 



trons. In this sub-section, we therefore zoom in on to 
the ti g bands and discuss the tight-binding hopping in- 
tegrals constructed out of the symmetrically orthonor- 
malized AMTOs for the truly minimal basis set of t^g 
orbitals. 

Although it would have been possible to compute the 
tight-binding Hamiltonian parameters in terms of three 



congruent orbitals d x > y i , d x ' z > , d ' z > , for the sake of com- 
parison with previous results, we preferred to work with 
the orbitals of eigen representation of symmetry lowered 
D?,d group, namely the a\ g and the two e g orbitals which 
transform according to a\ g and e g irreducible represen- 
tations. As mentioned in section II, the unitary trans- 
formation relating a\ g and e g with d x > v ', d x > z > , d x ' z > is 
given by, 

/ U a lg e*,l e£,2 \ 
d x > z > 1/V5 y/2/3 
dxV 1/V3 l/y/2 -l/s/6 

\ d y , z , 1/V3 -1/V2 -1/VeJ 

FIG.[5]shows the LDA band-structure of rhombohedral 
V2O3 plotted over an energy range of -1.5 eV below the 
Fermi level to 4 eV above the Fermi level, now projected 
on to cii g and e g degrees of freedom. e g (e g ) derived 
bands are also seen within the energy scale of the plot. 
We notice significant mixing between a l£( and e g charac- 
ters in the bands of interest spanning the energy range 
-1.0 to 1.7 eV, arising due to ai g -e g hopping processes 
between neighboring V sites. We also notice due to sym- 
metry reason e g orbitals acquire non-significant e g (e g ) 
character too. 

In the top left panel of FIG. [9] we show the downfolded 
band-structure obtained with truly minimal set consist- 
ing of ai g and e g . The bands, as they should be, are 
identical with the downfolded bands, shown in the second 
panel of FIG. El obtained using d x > y t, d x ' z > , d x < j min- 
imal basis set. The bottom left of the figure shows the 
dig and e g bands switching off the di g -e g hybridization. 
It is important to note that the width of the projected 
dig band in FIG.[5]is much more than that where the hy- 
bridization between d\ g and is neglected. Much of the 
d\g band width therefore comes from the hybridization 
with e g . This implies the crucial role of the d\ g and 
hybridization in proper description of the band-structure. 

The members of the truly minimal set are shown in 
FIG. [TUl the dig orbital which is oriented vertically and 
two more planar e g orbitals. Following the mixing be- 
tween e g (e g ) and e g as seen in the fatband plot of FIG. 
[51 we notice the significant presence of e g tails in the plots 
of e g orbitals, in addition to usual anti-bonding covalent 
character of the oxygen tails. The e g tails combine with 
the oxygen tails to produce sausage-like structures e.g. 
that at V sites located at 4 and 6 in the plot of , 2. 

Once we have defined the basis, in the following we 
compute the tight-binding hopping matrix elements be- 
tween the orthonormalized, truly minimal d\ g and e g 
AMTOs. As indicated in the Appendix, this is done 
by constructing H LDA (k.) in the Bloch k- representation, 
in the basis of symmetrically orthonormalized AMTOs, 
defined for di g and truly minimal set, for all 
k-points in the BZ and by subsequent Fourier transfor- 
mation H(k) — > H(r), for a given cluster with real-space 
range, r. Following CNR paper, for our real-space calcu- 
lation, we considered a cluster of fourteen V sites - the 
V sites belonging to such a cluster are marked in FIG. 
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O-p V-t2g V-eg V-s 




FIG. 7: Like in FIG. [5] but the solid lines in different panels now show the bands obtained with the downfolded 0-2p, V- 
t22 g , V-e 9 and V-4s LMTO basis in comparison to those obtained with full LMTO basis (shown in dotted lines). The I- and 
independent e„-s are shown on the right-hand side in each panel. 



[TJ The convention adopted for numbering the V atoms is 
same as that of CNR. All the V atoms which have at least 
one shared oxygen in the VOq octahedra with that of the 
central V atom (marked as 1 in FIG. [1]) are considered 
in this process. Distance- wise the short vertical VI- V4 
bond forms the nearest-neighbor (NN), while the three 
more or less planar bonds along the directions 2, 3, and 
3' form the 2nd neighbor shell with distances 6% larger 
than the V1-V4 bond-length of 2.697 A. The next shell 
of neighbors i.e. the 3rd neighbor shell is formed by the 
V atoms, 4, 8 and 8_ sitting at a distance 0.590 A farther 
than the second NN shell consisting of V2, V3, and V3 . 
The farthest shell of neighbors, 4-th neighbor shell in the 
cluster are formed by the V atoms, V5, V6, V6 , V5, V6, 
V<y whose distances from the central VI atoms differs by 
another 6% from that of V4, V8 and V8\ 

Before we actually attempt on numerically comput- 
ing the matrix elements of the tight-binding Hamilto- 
nian between the orthonormalized, downfolded a\ g and 
e g TVMTOs, a rough guess of the relative strength of var- 
ious hopping matrix elements may be obtained by exam- 
ining the overlap of TVMTOs placed at different V sites. 
In FIGs. fTTI and [T2l we consider two such representative 
cases where the truly minimal TVMTOs are placed along 
the vertical bond, 4, at sites VI and V4 (FIG. Q3]) and 
that along the horizontal bond, 2, at sites VI and V2 
(FIG. H"2"]) . The important feature to notice is the hop- 
ping paths via the oxygen and the e g tails in addition to 
the V-d — V-d hopping paths. While the approach taken 
by CNR, did considered the renormalization effect com- 
ing from oxygen degrees of freedom in some form, the 



effect due to the renormalization coming from e g 's was 
completed ignored, which has important consequences 
for hopping integrals e.g. connecting VI and V2 sites 
in FIG. H21 These additional hopping paths via the e g 
tails increase the importance of the hopping processes in 
the basal plane. Focusing on to the a\ g — a\ g vertical 
pair overlap, the biggest of all the hopping processes, we 
see the large, direct a\ g — a\ g hopping which is bonding 
{negative) in nature - the red lob at VI site overlaps with 
the red lob at V4 site. To understand the additional hop- 
ping contributions via the oxygen tails let us consider the 
renormalized a\ g orbitals at VI and V4 which consider- 
ing only the oxygen contributions can be written, from a 
simplistic point of view as: 

= |di)+A|pi) (2) 
|* 4 ) = \d 4 )+X\ P4 ) 

where A is the covalency mixing parameter between V-a\ g 
and O-p, and \p 4 ) are the wave- function of the shared 
O's between the VI Oq and V4C-6 octahedra having the 
same a\ g symmetry, \d\) and \di) are the bare a\ g orbitals 
at the VI and V4 sites. From this, we see that the overlap 
between the renormalized a\ g orbitals at sites VI and V4 
is given as , 

(*i,* 4 ) = {d 1 ,d A ) + \{(p 1 ,d i ) + (d u p 4 )) + X 2 (p llPi ) 

While both the d — d and p — p overlaps, {d\,di) and 
(pi,P4) are of bonding nature (the contribution (pi,P4) 
is small due to the presence of prefactor A 2 ), the sign 
of the correction terms, (pi,c?4) and (di,p4) depend on 
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PM 



FIG. 8: LDA band-structure of corundum V2O3 with the or- 
bital character projected on to a\ g and e g orbitals constructed 
out of three congruent t2 g orbitals, x y , x z and y z [see 
Eqn.(l)]. 



the V1-0-V4 angle, which for the real crystal turn out to 
be 82.3°. Careful investigation of FIG. [11] show these 
overlaps to be anti-bonding (positive)^. These anti- 
bonding hopping paths, therefore, oppose the bonding, 
direct a\ g — a\ g hopping, and thereby reduces the mag- 
nitude of the effective a\ g — a\ g hopping from the bare 
a ig—aig hopping. We will return to this point again while 
discussing the crucial sensitivity of this important hop- 
ping parameter on the intricate details of the geometrical 
structure. We further see that the direct vertical hopping 
processes between e g s are rather weak, which get weaker 
by the oxygen renormalization effect. 

Moving to FIG. [T^l for overlap along the horizontal 
bond direction, 2, we see a weak overlap between the 
ai g truly minimal TVMTOs while the e g like tails make 
the hopping between e g ,2 and e g ,2 nearly as strong as 
that between a\ g and e g , 2. The overlap e g , 1 — e g ,l is 
anti-bonding (positive) while e g , 2 — e g , 2 and a\ g — e g , 2 
overlaps are bonding (negative). 

In the first row of blocks in Table II we show all the 
hopping integrals between the central V atom (1) and 
the neighboring V atoms (for numbering sec FIG. 1) up 
to 4-th neighbor obtained by Fourier transform of the 
downfolded ti g Hamiltonian in a\ g — e g basis in symmet- 
rically orthogonalized representation. The on-site ener- 
gies are shown in Table I. We see that the hoppings be- 
yond the four predominant directions, 4, 2, 3, 3 also have 
non-negligible contributions. The tight-binding t2 g bands 
considering interactions till 2NN, 3NN and 4NN hoppings 







FIG. 9: Top panel: Downfolded band-structure of corundum- 
structured V2O3 in pure and Cr-doped phases obtained with 
the truly minimal basis set consisting of a\ g and e g ATVtTOs. 
The bands in the left are identical with the solid bands shown 
in the second panel of FIG. [5] and that of the bands obtained 
with the full NMTO basis in the energy range -1.5 below to 2 
eV above Ef=0. The members of the basis set are shown in 
FIG. 1101 Bottom panel: The a\ g (solid lines) and e g (dotted 
lines) bands in pure and Cr-doped V2O3 where the hybridiza- 
tion between a\ g and e g degrees of freedom has been switched 
off. Note the shrinkage of the total band-width compared to 
that in the left. 



are shown in FIG. [121 The fc-space band-structure con- 
sidering the infinite summation in the real-space Fourier 
series is also shown for comparison. As we see, hoppings 
till 4NN are essential to reproduce atleast the gross fea- 
tures of the band-structure. 

The symmetry properties of the corundum structure 
allows one to recast the important hopping integrals in 
the directions 4, 2, 3 and 3 in terms of reduced param- 
eters like (i, A, a, /3, a and r. The relationship of the 
various hopping integrals and the reduced parameters are 
shown in Table IIfi2. 

In CNR's original paper— as well as in the later paper 
by Di Matteo et. al^, the hopping integrals have been 
quoted in terms of these parameters. Di Matteo et. a/3 
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FIG. 10: Like in FIG. [6] but for the a lg and two e g iVMTOs, 
corresponding to downfolded bands in FIG. |U The ai 9 or- 
bital is oriented vertically while the e g orbitals are of more 
planar geometry. Due to symmetry reason, the e g orbitals 
bind strongly with e g degrees of freedom at other V sites 
while this binding is practically negligible for ai g . 



quoted the estimate obtained by CNR as well as the esti- 
mates obtained by TB fitting performed on LAPW band 
calculation of Mattheiss 3 -. These two estimates were 
found to be numerically not very different, although it 
was not clear whether the tight-binding fitting procedure 
included the correction due to the renormalization effect 
from the oxygen degrees of freedom or not. The esti- 
mate for the vertical pair hopping quoted by Di Matteo 
et. aZ. 13 appeared to be virtually same as the estimate 
of dda obtained by Mattheiss 36 , although one expects 
some difference due to integrated out oxygen degrees of 
freedom. For the sake of direct comparison and for the 
sake of future analysis, we have extracted these reduced 




(e|j(4)|H| e| 1 (i))=-o.06cv 




(e| 2 (4) \H\ e£ 2 (i)) = '0.06 e v 



FIG. 11: Overlap between V-ai 9 and V-e 9 downfolded 
TVMTOs, placed at the central VI site and the V4 site, forming 
vertical nearest-neighbor VI- V4 pair. The light(dark) shaded 
orbitals correspond to V1(V4). This gives the idea of impor- 
tant hopping paths. Note the hopping paths via the oxygen 
tails in addition to direct V-V paths. 



In Table IV we show the reduced parameters obtained 
by NWTO-downjolding technique in comparison to that 
of CNR and Di Matteo et. al. for the directions 4, 2, 3 
and 3. 



(a lg {i)\H\a lg (2))= 0.03 e v 




(e|j(i)|H] e£ | (2)) = +o.o8 ev 




(e| 2 0)|H| e| 2 (2)) = -o.2i ev 




FIG. 12: Same as FIG. [TT] but the TVMTOs are placed one 
at the central VI site and another at the neighboring V2 site 
along the x-axis. The light(dark) shaded orbitals correspond 
to VI (V2). Note the importance of e 9 tails in providing the 
hopping channels, in addition to oxygen mediated and direct 
V-V hopping channels. 
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1 .27 .00 .00 .00 .00 .00 .00 .00 .00 
| .30 .00 .00 .00 .00 .00 .00 .00 .00 
1 .28 .00 .00 .01 .01 -.04 -.04 .00 .00 

TABLE I: On-site matrix elements in the high-temperature 
paramagnetic metallic phase (undoped, ambient pressure 
V2O3), paramagnetic insulating ((Vo.962Cro. 038)203) and 
in the low-temperature antiferromagnetic insulating (mono- 
clinic) phase, between the m-orbital and the m'-orbital. Ex- 
cept for the orthonormalization, the orbitals are as defined in 
FIG. [10] 




FIG. 13: The tight-binding bands of corundum-structured V2O3 with various range of hopping interactions. In the first panel, 
hopping interactions include only the near neighbor interactions, interaction along the short vertical bond direction 4 and that 
along three basal bond directions, 2,3 and 3 . In the second panel interactions until third nearest neighbors (directions 4, 2, 3, 
3 , 4, 8 and 8_) and in the third panel interactions until forth nearest neighbors (directions 4, 2, 3, 3 , 4, 8, 8_, 5, 6, 6 , 5, 6 and 
6_) are included. The last panel shows the fully converged tight-binding band-structure involving all the hopping interactions 
ranging from nearest-neighbor to infinity. 



On examining the AMTO derived /1, A, a, 0, a and r 
parameters in comparison to that of CNR and Di Mat- 
teo et. al. we find that the parameters are quite different 
from their estimates. In particular we notice the signif- 
icant reduction of the vertical pair a% g — a\ g hopping, 
given by the parameter p and the increased importance 
of the hoppings in the basal plane. In order to inves- 
tigate the probable reasons for such discrepancy, one of 
the prime candidate on first glance appears to be the so- 
phisticated treatment of AMTO-downfolding over that 
of CNR, where the effective orbitals were constructed 



following Anderson's super-exchange idea, the covalency 
V-tf— O-p mixing parameter, A and the charge transfer 
gap E^d — Ei v were extracted from nuclear magnetic res- 
onance and photo-emission experimental measurements 
which were known only to certain accuracy, and finally 
the matrix elements were computed in terms of second 
order perturbation theory in A. Nevertheless, in spite of 
all the above mentioned approximations - some of which 
are crude- it turns out that the structural information 
plays a even more crucial role. 

We consider in the following the specific case of verti- 
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TABLE III: Hopping parameters along the directions 4,2, 
3 and 3 in terms of corundum-symmetry-adopted reduced 
parameters fi, p, a, f3, a and r. 
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TABLE IV: Comparison of TB parameters of corundum- 
structured V2O3 (in terms reduced parameters fi, p, a, f3, 
a and r) obtained by different procedures for V-V hopping 
along the four near directions, 4, 2, 3 and 3. 

cal a\ g — aig hopping which doesn't have the additional 
complexity of hopping via the e g tails, another crucial in- 
gredient not taken into account in CNR's study. In CNR 
paper, the direct d — d hoppings were obtained from a 
linear combination of atomic orbital (LCAO) kind of ap- 
proach by Ashkenazi et. a4£ which assumes the correct 
geometry, while the trigonal distortion was assumed to 
negligible (set to zero)^. Such an approximation is found 
to have a deeper implication in terms of the quantitative 
estimates of the effective V — V hopping. In order to 
have an understanding of the delicate effect of the ge- 
ometry, we carried out calculations on crystal structures 
with varying trigonal distortions. Crystal structures with 
varying amount of trigonal distortions are generated by 
linear interpolation of the internal parameters associated 
with V and O atoms between that of the real crystal and 
that of the ideal hexagonal arrangement: 

z v = 0.3333(1 -c) + 0.34630c 
x = 0.3333(1 -c) + 0.31164c 

Putting c = 0(1) gives the ideal(real) structure. In- 
creasing c increases the trigonal distortion which is the 
difference between two sets of V-0 distances in VC>6 octa- 
hedra, with three short and three long V-0 bond-lengths. 
Changing the parameter c, however also changes the VI- 
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Trigonal distortion ( A°) 
0.00 0.02 0.04 0.06 0.08 



Real siructure 



Ideal stivcture 




V-V dist. (A u ) 

FIG. 14: The influence of the trigonal distortion and VI- V4 
distance on the vertical a\ g — ai g pair hopping magnitude. 
Structures with varying amount of distortions are generated 
by linear interpolation of the internal parameters associated 
with V and O positions, z v = 0.3333(1 — c) + 0.34630c; zo ~ 
0.3333(1 — c)+0. 311640c. This leads to simultaneous variation 
of the VI- V4 distance and the trigonal distortion. The dotted 
and the solid lines give the estimate of the direct oi s — oi s 
hopping and that of the effective ai g — oi s hopping in truly 
minimal a% g — e g 7VMTO basis. The open and solid circles 
are the estimates of direct and effective ai g — ai g hopping 
respectively, as obtained by CNR 4 . The VI- V4 distance in 
CNR's calculation were considered to be same as that in real 
structure while the amount of trigonal distortion was set to 
zero. 



V4 distance which effects the bare or direct a± g — a± g 
hopping. We have carried out calculations for c=0, 1/3, 
2/3 an 1. The calculations for the hopping matrix el- 
ements are carried out for the truly minimal set of a\ g 
and e g iVMTOs as well as for the set where the O-p's 
are kept active in addition to V-d's. This has been done 
to bring out the renormalization effect coming from oxy- 
gen degrees of freedom. The results are shown in FIG. 
[HI We see that the pd contribution given by the differ- 
ence of direct a\ g — a± g and the renormalized a\ g — a\ g 
hopping, almost vanishes for the ideal structure and in- 
creases monotonically as the trigonal distortion increases 
towards to the value obtained in the real structure. The 
pd contribution, apart from the case of ideal structure, is 
anti-bonding whereas the direct (bare) dd contribution is 
bonding. On top of this, comes the even stronger trend 
that the bonding dd interaction decreases with the VI- 
V4 bond distance, d, as d -3,3 . As a result, a\ g — a\ g 
hopping integral depends strongly on the distortion. As 




-U«cV 



FIG. 15: Comparison between Vl-ai 9 and V4-ai 9 overlaps for 
the real and ideal structures. The light(dark) shaded orbitals 
correspond to V1(V2). The tails of the orbital at VI site has 
been omitted for clarity. 



already mentioned, in CNR's calculation, though the 
bare dd hopping was obtained following LCAO type of 
approach^ on a correct V-V geometry with correct VI- 
V4 bond distance, the trigonal distortion of the VOe oc- 
tahcdra was assumed to be negligible with all V-0 bond 
lengths to be equal. This resulted into a direct dd hop- 
ping of -0.66 eV in good agreement with our correspond- 
ing estimate of -0.72 eV, while the pd contribution gave 
rise to a small renormalization of -0.06 eV. We note that 
the pd contribution is small and of bonding (negative) 
type as we obtained in our calculation with ideal struc- 
ture. This in turn proves the extreme sensitivity of the 
hopping parameters on the correct geometry of the sys- 
tem. In FIG. [T5]we show the comparison between over- 
lap of Vl-ai 9 and V4-ai y ATVITOs for the real and ideal 
structures. Studying the figure, we notice that due to the 
about 14 % reduction of V1-V4 distance in case of ideal 
structure, the direct, bonding type ai g — a± g overlap is 
much stronger than compared to that in real structure 
which gives rise to the strong slope observed in FIG. Q~J] 
Focusing on the pd contribution, for the ideal structure, 
the p— type tail from V4 MvITO at 05 site (referred as 
Pi in Eqn.2) passes almost through the node of the a\ g 
orbital at VI site, as shown with solid line, and gives 
rise to a negligible overlap between d\ and p^ (c.f. Eqn. 
[3]). For real structure, on the other hand, the change of 
V1-05-V4 angle causes the pi tail having finite overlap 
with a\ g orbital at VI site which turn out to be positive 
(antibonding) in sign. 
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D. Vertical pair model and validity of molecular 
orbital states 

In this sub-section, we discuss the consequences of the 
new set of JVMTO derived parameter values in the con- 
text of validity of vertical pair model and the molecular 
orbital formation. 

Di Matteo et. alJ£ analyzed in detail the parameter 
space of the V2O3 problem in context of all possible or- 
bital and magnetic ground states configuration of the ef- 
fective many-body Hamiltonian by using variational pro- 
cedure. The orbital wave- functions of the ferromagnetic 
state of the vertical pair of V atoms at sites a and b were 
postulated to be given by: 

|V±) afc - 1/V2(| ± l) o |0) 6 + I ± 1) 6 |0)„) (3) 

where |0) = |e* lej,2), | - 1) = \a lg e*,l) and |1> = 
\ai g eg, 2) are three two-electron states constructed out of 
three one-electron states, e£ , 1, e^, 2 and a\ g . 

The correlation energy of such a state was defined as 
the difference between the ground state energy and its 
Hartree-Fock approximation. Taking into account the 

hopping integrals in terms of reduced parameters, the 

( — 

ground state energy of such a pair is given by — P J , 
which involves back and forth virtual hoppings out of 

e g state [—jfz-j], that out of a\ g state [— jfzj] an d the 
correlated hopping between V atoms at a and b sites, 
where they exchange electrons in a lg and e" g states simul- 
taneously. U and J are Coulomb and exchange integrals 
respectively. The latter mechanism which arises due to 
entangled nature of the form of the wave-function in (3) is 
absent in its Hartree-Fock approximation and gives what 
Di Matteo et. al. termed as molecular correlation energy. 
Two different regimes of solutions were defined depend- 
ing on the relative magnitudes of the correlation energy 
of the ferromagnetic state of the vertical pair, namely, the 
molecular correlation energy and the in-plane exchange 
energy (governed by the hopping processes in the basal 
plane) : 

(a) If the molecular correlation energy is larger than 
the in-plane exchange energy then the whole crystal con- 
sists of some molecular units and the variational wave- 
function should be constructed in terms of molecular 
states given by (3). 

(b) If the in-plane exchange energy is larger than the 
molecular correlation energy which tend to break the 
stability of the correlated molecular states, one needs to 
construct the variational wave- function in terms of single 
site atomic states. 

Di Matteo et. al. 13 approximated the in- plane ex- 
change energy by pa ^jfzrj ' [considering the hoppings 

along three 2NN bonds in the direction of 2, 3 and 3 
and neglecting the hoppings involving reduced parame- 
ters (3 and a] 

Considering the numerical values of /i, p, a and r, as 
given by CNR and Di Matteo et. al. one gets: 



2pp/{a 2 +r 2 ) Pa 4- 3.5 

while using the parameters obtained by iVMTO- 
downfolding one gets: 

2pp/{a 2 +t 2 ) re .83 

[this ratio gets even smaller taking into account the 
hoppings related to /3 parameter which is almost as large 
as r.] 

Therefore, while the TB parameters given by CNR and 
Di Matteo et. al. favor the formation of stable molec- 
ular orbital states, the parameters provided by iVMTO- 
downfolding procedure clearly do not favor it. We there- 
fore believe that it is very much needed to repeat the 
calculations using the new set of hopping parameters to 
shed light on the long standing puzzles in V2O3. This 
issue has been recently taken up in Ref. 41 . 

E. Cr-doped V2O3 ((Vo^Cro.cas^Os): the 
low-energy, tight-binding Hamiltonian 

The right, top panel of FIG. 9 shows the LDA band 
structure of 3.8% Cr doped V2O3 in the paramagnetic 
insulating phase. Upon comparison with the band struc- 
ture of the undoped V2O3 in the paramagnetic metal- 
lic phase, as presented in the left, top panel of FIG. 
9, one finds that the t2 S bandwidth is reduced in the 
doped compound to about 2.25 eV from about 2.5 eV 
in case of undoped compound. The right, bottom panel 
of FIG. 9 shows the ai s and e g bands switching off the 
aig-e^ hybridization. The lattice expansion upon Cr dop- 
ing, causes bottom of the ai ff band to move up and the 
top of e g band to move down, giving rise to max{e eg } 
- min{e a i g } pa 1.69 eV, in comparison to max{e eg } - 
min{e Q i g } pa 2.02 eV in case of undoped compound. 

The middle rows of Table I and II, lists the onsite en- 
ergies and hopping matrix elements corresponding to Cr 
doped V2O3. The crystal field splitting is found to in- 
crease by 0.03 eV compared to undoped case. The magni- 
tude of the dominant V1-V4 hopping is found to decrease 
from 0.51 eV in the undoped case to 0.43 eV in the doped 
case. These changes in the one-electron parameters were 
found to be significant to drive the metal-insulator tran- 
sition as explained in Ref 28 . 

F. Monoclinic V2O3: the low-energy, tight-binding 
Hamiltonian 

Finally, we thought it will be worthwhile to study 
the influence of the monoclinic distortion in the low- 
temperature crystal structure on the hopping integrals. 
It is of interest to know how much the crystal struc- 
ture change effects the hopping matrix elements. We ap- 
plied the same AMTO-downfolding machinery, described 
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in great detail in previous sections for V2O3 in corun- 
dum structure. The LDA self-consistent potentials are 
generated by TB-LMTO-ASA calculation with poten- 
tial sphere overlap less than 18 % and empty sphere 
overlap less than 22 %. The truly minimal a\ g and 
e g basis sets are defined within the framework of the 
iVMTO-downfolding technique for the monoclinic struc- 
ture. For the sake of comparison, we retained the 
corundum-symmetry-adopted a\ g and e g basis also in the 
monoclinic structure. However, the further lowering of 
the symmetry in the monoclinic phase introduces mixing 
between a\ g and e g orbitals at the same site, which is 
reflected as crystal field terms in the on-site block of the 
real-space Hamiltonian. The monoclinic distortion makes 
the various near neighbor distances unequal^ compared 
to that in corundum phase. While the vertical V1-V4 
bond expands by 1.8 %, the horizontal V1-V2 bond ex- 
pands by about 4 % making the three, basal near neigh- 
bor bond distances along 2, 3 and 3 unequal. It also 
makes the bond distance along 4 different from those 
along 8 and ST, the bond distances along 5 and 5 different 
from that of 6 and 6 and, 6 and 6_. The tight-binding 
hopping integrals and hopping elements computed as ele- 
ments of the orthonormalizcd a\ g and ATVITO Hamil- 
tonian are quoted in Table I and II. Focusing on the 
o-ig ~ a ig vertical pair hopping, we find that the value is 
further decreased to -.44 eV compared to the value of -.51 
eV in the corundum structure. This reduction is primar- 
ily driven by the 1.8 % increase in the VI- V4 bond length 
and slight tilting of the V1-V4 bond which decreases the 
magnitude of the bare dd hopping from -.72 eV in the 
corundum structure to -.64 eV in the monoclinic struc- 
ture. The pd contribution due to integrated out oxygen 
tails turned out to be +0.20 eV which can be compared 
with the value +.21 eV, that in corundum structure. The 
V1-0-V4 angle remains essentially unaltered between the 
corundum and monoclinic structure. 

Comparing the hoppings in other directions, the reduc- 
tion is maximum for the horizontal V1-V2 bond which 
expands by 4% over the value in corundum structure. 
The other two near neighbor bonds in the basal plane, 3 
and 3 on the other hand contracts (the bond 3 by 0.7 
% and 3 by .2 %) which is reflected in the changes in 
hopping integrals. Similar distance dependent increase 
or decrease can be observed for farther ranged hoppings. 

In brief, though the low-temperature structural change 
induces changes in the hopping parameters, these 
changes are not drastic. Till date, a direct, experimen- 
tal evidence of orbital ordering is lacking and issue of 
orbital ordering, its existence and type in V2O3 still re- 
mains highly controversial 42 . Nevertheless, the role of 
orbital degrees of freedom in stabilizing 5 the magnetic 
structure with broken trigonal symmetry of the corun- 
dum lattice remains to be plausible idea. In that case, 
the monoclinic distortion may possibly be the reflection 
of the peculiar spin and orbital ordering rather than the 
cause for it. 



IV. SUMMARY AND OUTLOOK 

To summarize, employing the NMTO-downfolding 
technique and Wannier function representation of the 
Hamiltonian, we have derived in a first-principles man- 
ner, the effective V-V hopping interactions correspond- 
ing to the low-energy, t 2g bands of V2O3. Our results 
show, contrary to popular believe, for modeling of V2O3, 
inter-pair V-V hoppings are equally important as V-V 
intra-pair hoppings. The significant changes in hopping 
parameters compared to CNR parameters occur primar- 
ily due to the neglect of trigonal distortion in the pre- 
vious study and due to the hopping processes via the 
integrated out e g tails in addition to that via oxygen like 
tails, a fact not considered before. This calls the need for 
revisiting the many-body calculations, which start with 
the assumption of vertical pairs as the building blocks. 

The Wannier functions corresponding to the low- 
energy, t2 9 bands derived in this paper will serve as the 
basis to define low-energy, multi-orbital Hubbard Hamil- 
tonian for the LDA+DMFT calculations for V2O3, which 
rely on the choice of flexible, atom-centered, localized ba- 
sis sets. Such calculations have been already carried out. 
For details please see Re@$. Considering the rather de- 
localized nature of the real-space Hamiltonian of V2O3, 
it is quite natural to expect improvements on going be- 
yond the single-site approximation of DMFT and tak- 
ing into account the cluster effect. For such study, it 
is crucial to decide on a minimal cluster which has the 
dominating effect and the vertical pair has been often dis- 
cussed as a natural choice. However, in view of iVMTO- 
downfolding derived parameters and the breakdown of 
correlated molecular orbital like states, the most tempt- 
ing choice of the vertical pair as the cluster seems to 
be hardly satisfactory. The cluster LDA+DMFT calcu- 
lations with VI- V4 pair show no qualitative difference 
with single site DMFT results for corundum PI phase 43 . 
The minimal cluster should include in addition to verti- 
cal pairs VI and V4, the near-neighbor V atoms in the 
basal plane, V2, V3 and V3 , which though is a computer 
expensive DMFT job to carry out. 
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APPENDIX A: 7VMTO METHOD 

In the following, we describe the AMTO method which 
provides a tool for direct generation of localized Wan- 
nier functions. The downfolding procedure is also im- 
plemented in this framework to construct truly minimal 
basis sets which pick out selected bands. 

In the ATVITO method, a basis set of localized orbitals 
is constructed from the exact scattering solutions for a 
superpositions, J2b, Vr ( Tr )^ °^ short-ranged, spherically- 
symmetric potential wells - a so-called MT approxima- 
tion to the potential. This is done by first numeri- 
cally solving the radial Schrodinger's equations, to find 
<Pri ( £ n, tr) Yi m (tr), the partial waves, for all angular 
momenta, I, with non- vanishing phase-shifts, for all po- 
tential wells, R, and for a chosen set of energies spanning 
the region of interest, e„ = eo, £n- 

- [ripm (e, r)]" =[e- v R (r) -1(1+1) /r 2 ] np m (e, r) 

The partial-wave channels, Rim are partitioned into 
active and passive channels. The active channels are 
those for which one chooses to have orbitals in the ba- 
sis set, i.e. they are the chosen one-electron degrees of 
freedom. The passive channels are said to be downfolded. 

For each active channel, Rim, a so-called kinked par- 
tial wave (KPW), (frgjm (e„,r) , is constructed. A kinked 
partial wave is basically a partial wave with a tail joined 
continuously to it with a kink at a central, so-called hard 
sphere of radius clr. The tail of the kinked partial wave 
is a screened spherical wave, ipRim ( £ j r ) i which is essen- 
tially the solution with energy s of the wave equation in 
the interstitial between the hard spheres, 

-A^> (e,r) = eip (e,r) 

with the boundary condition that, independent of the en- 
ergy, ipRim (e, r) go to Yj m (r R ) at the central hard sphere, 
and to zero (with a kink) at all other hard spheres at the 
neighboring sites. It is this latter confinement, which 
makes the screened spherical waves and the KPWs lo- 
calized when the energy is not too high. The default 
value of the hard-sphere radii, or, is 90 % of the ap- 
propriate covalent, atomic, or ionic radius. The above- 
mentioned boundary condition only applies to the active 
components of the spherical-harmonics expansions of the 
screened spherical wave on the hard spheres. For the re- 
maining downfolded or passive components the screened 
spherical wave equals the corresponding partial-wave so- 
lution of Schrodinger's equation throughout the MT- 
sphere, i.e. it has the proper phase shift. 

If one can now form a linear combination of such kinked 
partial waves with the property that all kinks cancel, one 
finds a solution of Schrodingers equation with energy e„. 
In fact, this kink-cancellation condition leads to the clas- 
sical method of Korringa, Kohn and Rostoker 44 (KKR), 
but in a general -so-called screened- representation and 
valid for overlapping MT potentials to leading order in 
the potential overlap. The screened KKR equations are a 



set of energy-dependent, homogeneous linear equations, 
with a matrix, Ksr~ ar- (s) , whose rows and columns 

' Him, Kim v > 5 

are labeled by the active channels. In the AMTO 
method, we don't solve this set of secular equations, but 
proceed via construction of energy- independent MTO ba- 
sis sets which span the solutions &i (r) with energies 
of Schrodinger's equation to within errors proportional 
to (e l - eo) (e» - ei) .. (e 4 - ejv) , where e , ei, e N is the 
chosen energy mesh with N+l points defined already. 
Such an energy-independent set of Nth-order MTOs is 
called an AMTO set. 

The members of the ATVITO basis set for the energy 
mesh eo, ejv are superpositions, 

N 

X W (r) = Y Y <j> nT ^{e n ,r)L {N l __ (Al) 

A-Rlm^ I L^i L^i ^Rlm V ™> I n Rlm,Rlm y ' 

n=0 Rim 

of the kinked partial waves, (f'mm ( e > r ) > a * the N + l 
points (labeled by n) of the energy mesh. Expression (1) 
is the energy-quantized form of Lagrange interpolation, 

N N 

X W( £ )«^0(e„)^)( £ ), J] ^-^ 

of a function of energy, <fi (e) , by an TV th-degree polyno- 
mial, x^ N ^ i e ) '■ The Ath-degree polynomial, (e) , is 
substituted by a matrix with elements, L^ N )- — , the 

J ' nRlm,Rlm ' 

function of energy, <j> (e) , by a Hilbert space with axes, 
4>Rlm ( e > r ) ' an d ^e interpolating polynomial, x ( £ ) , 
by a Hilbert space with axes, X j$ m ( r ) ■ 

Note that the size of the AMTO basis is given by the 
number of active channels and is independent of the num- 
ber, N + 1, of energy points. The energy-selective and 
localized nature of AMTO basis makes the AMTO set 
flexible and may be chosen as truly minimal, that is, to 
span selected bands with as many (few) basis functions as 
there are bands. If those bands are isolated, the AMTO 
set spans the Hilbert space of the Wannier functions and 
the orthonormalized AMTOs are the Wannier functions. 
But even if the bands of interest overlap other bands, it 
may be possible to pick out those few bands and their 
corresponding Wannier-like functions with the AMTO 
method. The AMTO method can thus be used for di- 
rect generation of Wannier or Wannier-like functions. 

The Lagrange coefficients, as well as the Hamil- 

tonian and overlap matrices in the AMTO basis are ex- 
pressed solely in terms of the KKR resolvent, K (e) 1 , 
and its first energy derivative, K (e) , evaluated at the 
energy mesh, e = eo, cat. Variational estimates of the 
one-electron energies, Si, may be obtained from the gen- 
eralized eigenvalue problem, 




with 

H = -A + ^v r {\t-R\), 

R 
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or as the eigenvalues of the one-electron Hamiltonian ma- 
trix, 



H 



LDA 



XN)± 



\H\X 



(N)J 



in the basis of symmetrically orthonormalized MvITOs: 



(N)± 



X y 'J \T ' \ X y i 



In the present paper, the orbitals shown are TVMTOs 
before orthonormalization because they are (slightly) 
more localized than the orthonormalized ones. The hop- 
ping integrals and on-site elements given in the tables are 
of course matrix elements of the orthonormalized Hamil- 
tonian. 

For crystals, all calculations except the generation of 
the screened structure matrix are performed in the Bloch 
k- representation 



V 



T 



T)exp{27rik- (R + T)} 



where T labels the L(— > oo) lattice translations and R 
the active sites in the primitive cell. In order to ob- 
tain the orbitals and the Hamiltonian in configuration 
space, Fourier-transformation over the Brillouin zone is 
performed. 

It is worth-mentioning here that this construction of 
a minimal ATVITO basis set is different from standard 
Lowdin downfolding. The latter partitions a given, large 
(say orthonormal) basis into active (A) and passive (P) 
subsets, then finds the downfolded Hamiltonian matrix 
as: 

(A(e)\H\A(e)) = (A \H\ A) 

- (A \H\ P) (P \H - e\ Py 1 (P \H\ A) 

and finally removes the e-dependence of the downfolded 
basis by linearizing (P\H — e\P)~ 1 and treating the 
term linear in e as an overlap matrix. Obviously, since 
the AMTO set is exact at N+ 1 energy points, it is more 
accurate. 

Our present MvITO code is however not yet self- 
consistent, so we used the current Stuttgart tight-binding 
version of the linear- muffin-tin-orbital (TB-LMTO) 34 
code within the atomic sphere approximation (ASA) to 
generate the LDA potentials. Despite this shape approx- 
imation for the potential, the iVMTO bands used in the 
present paper are more accurate than LMTO bands, first 
of all because the ./VMTOs do not use the zero-energy ap- 
proximation in the interstitial region and, secondly, be- 
cause we use N > 1. 



APPENDIX B: COMPUTATIONAL DETAILS 

As mentioned in Appendix A, our present iVMTO 
code is not self-consistent, we therefore used the current 



Stuttgart TB-LMTO-ASA code to generate the LDA po- 
tentials. Such a potential in the atomic-spheres approxi- 
mation is an overlapping MT-potential, like the one han- 
dled by the iVMTO method, but with the relative over- 
laps, 



IR-R'I 



1, 



(Bl) 



limited to about 20%. This limitation comes from 
the LMTO-ASA+cc method, which solves Schrodinger's 
equation by treating the overlap as a perturbation (the 
so called combined-correction term, cc) and uses screened 
spherical waves of zero kinetic energy in the s-interstitial. 
Poisson's equation is solved for the output charge den- 
sity, spherically symmetrized inside the same atomic s- 
spheres. 

We now specify our computational set-up. The radii 
of the potential spheres, s#, were dictated by our use 
of the LMTO-ASA method to generate the LDA poten- 
tials. In order to limit the overlaps defined by equation 
(IBip . interstital -or empty- spheres (E) were inserted 
in the non-cubic structures. Table [V] gives the radii of 
the potential spheres. As a result, the overlap between 
atomic spheres was <16%, between atomic and empty 
spheres <18%, and between empty spheres <20%. We 
used the guidance given by the current version of the 
code in choosing the potential spheres appropriately. 

TABLE V: Radii sr of potential spheres in Bohr atomic units. 
V Ol 02 E El E2 E3 E4 E5 



PM 2.46 1. 
PI 2.46 L 
AFI 2.46 L 



2.36 1.80 
2.36 1.80 

2.36 2.36 1.80 1.80 1.80 1.80 



TABLE VI: LMTO basis sets used in the self-consistent cal- 
culation of LDA potential. (I) means that the i-partial waves 
were downfolded within in the LMTO-ASA+cc. 



V Ol 



02 



E 



El 



E2 E3 E4 E5 



PM spd sp(d) s(pd) s(p) 

PI spd sp(d) s(pd) s(p) 

AFI spd sp(d) sp(d) s(pd) s(pd) s(p) s(p) s(p) s(p) 



The self-consistent valence-electron densities were cal- 
culated with the LMTO bases listed in Table ED We 
found it is important to downfold the oxygen d partial 
waves, rather than to neglect them (i.e. to approxi- 
mate them by spherical Bessel functions when solving 
Schrodingers equation, and to neglect them in the charge 
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density). Since the LMTO calculations were used to pro- 
duced the self-consistent charge densities, the energies, 
em , for the linear 4>m , <j)m expansions were chosen at the 
centers of gravity of the occupied parts of the respective 
DOS ^-projections. 

Finally, in the ATVITO calculations, the hard-sphere 
radii, an, for the active channels were chosen as 0.7sr. 
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with that obtained by LAPW calculation of Matteheiss. 
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blue pa lob overlap with red d\ lob). 
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followed by translation, if necessary. This, for example, 
introduces a negative sign for the e^,l wavefunction at 
site V4 compared to that at site VI. In CNR paper, they 
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